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ABSTRACT 

We consider the general conditions of quark droplets formation in high density 
neutron matter. The growth of the quark bubble (assumed to contain a suffi- 
ciently large number of particles) can be described by means of a Fokker-Planck 
equation. The dynamics of the nucleation essentially depends on the physical 
properties of the medium it takes place. The conditions for quark bubble forma- 
tion are analyzed within the frameworks of both dissipative and non-dissipative 
(with zero bulk and shear viscosity coefficients) approaches. The conversion time 
of the neutron star to a quark star is obtained as a function of the equation of 
state of the neutron matter and of the microscopic parameters of the quark nu- 
clei. As an application of the obtained formalism we analyze the first order phase 
transition from neutron matter to quark matter in rapidly rotating neutron stars 
cores, triggered by the gravitational energy released during the spinning down of 
the neutron star. The endothermic conversion process, via gravitational energy 
absorption, could take place, in a very short time interval, of the order of few 
tens seconds, in a class of dense compact objects, with very high magnetic fields, 
called magnet ars. 

Subject headings: dense matter - pulsar : general - stars : interiors - stars : 
neutron 



1. Introduction 

Since Witten (1984), following early proposals by Itoh (1970) and Bodmer (1971), sug- 
gested that strange quark matter, consisting of u-, d- and s-quarks is energetically the most 
favorable state of the matter, the problem of existence of strange quark stars has been 
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intensively investigated in the physical and astrophysical literature. The possibility that 
some compact objects could be strange stars remains an interesting and intriguing, but still 
open question. Witten (1984) also proposed two ways of formation of strange matter: the 
quark-hadron phase transition in the early universe and conversion of neutron stars into 
strange ones at ultrahigh densities. In the theories of strong interaction quark bag models 
suppose that breaking of physical vacuum takes place inside hadrons. As a result vacuum 
energy densities inside and outside a hadron become essentially different, and the vacuum 
pressure on the bag wall equilibrates the pressure of quarks, thus stabilizing the system. If 
the hypothesis of the quark matter is true, then some of neutrons stars could actually be 
strange stars, built entirely of strange matter (Alcock, Farhi & Ohnto 1986; Haensel, Zdunik 
& Schaeffer 1986)). However, there are general arguments against the existence of strange 
stars (Caldwell & Friedman 1991). For an extensive review of strange star properties see 
Cheng, Dai & Lu (1998). 

Several mechanisms have been proposed for the formation of quark stars. Quark stars 
are expected to form during the collapse of the core of a massive star after the supernova 
explosion, as a result of a first or second order phase transition, resulting in deconfined quark 
matter (Dai, Peng & Lu 1995). The proto-neutron star core, or the neutron star core, is a 
favorable environment for the conversion of ordinary matter to strange quark matter (Cheng 
et al. 1998). Another possibility is that some neutron stars in low-mass X-ray binaries can 
accrete sufficient mass to undergo a phase transition to become strange stars (Cheng & 
Dai 1996). This mechanism has also been proposed as a source of radiation emission for 
cosmological 7-ray bursts (Cheng & Dai 1998). Some basic properties of strange stars hke 
mass, radius, cooling, collapse and surface radiation have been also studied (Cheng & Harko 
2000, 2003). 

The physical mechanisms of the transition from neutron matter to quark matter in 
an astrophysical background have been studied within several models. The first is due to 
Olinto (1987), who used a non-relativistic diffusion model. As such, this is a slow combustion 
model, with the burning front propagating at a speed of approximately 10 m/scc. This is 
determined primarily by the rate at which one of the down quarks inside the neutrons is 
converted, through a weak decay, to a strange quark: d + u ^ s + u. The second method 
of describing the conversion process was first suggested by Horvath & Benvenuto (1988), 
and analyzed in detail by Lugones, Benvenuto & Vucetich (1994) and Lugones & Benvenuto 
(1995), who modelled the conversion as a detonation. In this case the conversion rate is 
several orders of magnitude faster than that predicted by the slow combustion model. This 
model is based on the relativistic shock waves and combustion theory. But regardless of the 
way in which the transformation occurs, an initial seed of quark matter is needed to start 
the process. 
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Real neutron stars have two conserved charges- electric and baryonic. Therefore, a neu- 
tron star has more than one independent component and in this sense it is a "complex" 
system (Glendenning 1992). The characteristics of a first-order phase transition, like the 
dcconfinement transition, are very different in the two cases. In the "complex" system the 
conserved charges can be shared by the two phases in equilibrium in different concentrations 
in each phase. The mixed phase, formed from hadrons and quarks cannot exist in a simple 
body in the presence of the gravity because the pressure in that phase is a constant (Glen- 
denning 2000). This causes a discontinuity in the density distribution in the star occurring 
at the radius where Gibbs criteria are satisfied. The isospin symmetry energy in neutron rich 
matter will exploit the degree of freedom of readjusting the charges between hadronic and 
quark phases in equilibrium so as to reduce the symmetry energy to the extent consistent 
with charge conservation. Regions of hadronic matter will have a net positive charge neu- 
tralized by a net negative charge on the quark matter regions (Glendenning 2000). Coulomb 
repulsion will prevent the regions of like charge to grow too large, and the surface energy 
will act in the opposite sense in preventing them to become too small. Consequently, the 
mixed phase will form a Coulomb lattice so as to minimize the sum of Coulomb and surface 
interface energy at each proportion of phases. As the quark phase becomes more abundant, 
the droplets merge to form rods, rods merge to form slabs etc. (Glendenning 2000). Hence 
the actual geometric phase of quark nuggets in a neutron star evolves as a function of the 
pressure and of the surface energy between dense nuclear matter and quark matter. 

In order to describe the process of formation and evolution of quark nuclei in the neutron 
stars one must use nucleation theory. The goal of nucleation theory is to compute the 
probability that a bubble or droplet of the A phase appears in a system initially in the 

B phase near the critical temperature (Landau & Lifshitz 1980). Homogeneous nucleation 
theory applies when the system is pure. Nucleation theory is applicable for first-order phase 
transitions when the matter is not dramatically supercooled or superheated. If substantial 
supercooling or superheating is present, or if the phase transition is second order, then the 
relevant dynamics is spinodal decomposition (Shukla & Mohanty 2001). 

The nucleation of strange quark matter inside hot, dense nuclear matter was investi- 
gated, with the use of Zel'dovich's kinetic theory of nucleation, by Horvath, Benvenuto & 
Vucetich (1992). By assuming that the newly formed strange quark matter bubble can be 
described by a simple bag model containing Ng ultrarelativistic quarks, and by assuming 
that the time evolution of the radius of the bubble is described by an equation of the form 
dr/dt — {r — Rc)/T^, where Rc is the critical radius of the bubble and the weak interaction 
time scale, the nucleation rate ^ is given by 

e = 2.2 X 10-2 {T/af' N'J%' exp (-3.1iVVV/T) , (1) 
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where Nqc < 300 is the critical quark number, a is the surface tension of the bubble and 

Tin and Uq are the particle number densities in the neutron and quark phase, respectively 
(Horvath, Benvenuto & Vucetich 1992). In this way a lower bound for the temperature of 
the nucleation, T > 2.1MeV, can be obtained. The effects of the curvature energy term 
on thermal strange quark matter nucleation have been considered by Horvath (1994), who 
derived, within the same approach, the following expression for the nucleation rate: 



All the effects of the curvature enter only through the value of Re- As a general conclu- 
sion one obtains the result that even the curvature term acts against strange quark matter 
nucleation, the physical temperature of a just born proto-neutron star is, in any model, more 
than enough to drive an efficient boiling of the neutron material (the observations of the neu- 
trino flux from SN1987A are consistent with an effective temperature of T ~ 4 MeV). The 
possibility of stable strange quark matter (both bulk and quasi-bulk) at finite temperature, 
and some of its properties have been investigated, within the framework of the dynamical 
density-dependent quark mass model of confinement, by Chakrabarty (1993). The behavior 
of the surface tension and the stability of quark droplets at T 7^ has also been discussed, 
with reference to strangelet formation in ultra-relativistic heavy-ion collisions. 

A different approach to the problem of thermal nucleation was taken by Olesen & 
Madsen (1994). They used the assumption (standard in the theory of bubble nucleation in 
first-order phase transitions) that bubbles form at a rate given by it! exp (— Wc/T), 

with Wc the minimum work required to form a bubble with radius Tc- As a main result 
it follows that if the bag constant lies in the interval where three-flavor but not two-flavor 
quark matter is stable at zero pressure and temperature (145 MeV< B^^^ < 163 MeV), then 
all or parts of a neutron star will be converted into strange matter, during the first second of 
its existence. For bag constant above the stability interval, a partial transformation is still 
possible (Olesen & Madsen 1994). 

Heiselberg (1995) calculated the rate of formation of quark matter droplets in neutron 
stars from a combination of bubble formation rates in cold degenerate and high tempera- 
ture matter, taking into account nuclear matter calculations of the viscosity and thermal 
conductivity. The droplet formation rate is that of Langer & Turski (1973), given by / = 
{k/2'K) Q.Q exp {—Wc/T), with fig = (2/3\/3) (cr/T)^^^ {fc/^q}^ the statistical prefactor, which 
measures the phase-space volume of the saddle point around Rc that the droplet has to pass, 
in its way to the lower energy state, is the quark correlation length. The dynamical prefac- 
tor determines the droplet growth rate and is given by A; = {2a / Aw'^r^) {XT + 2 (477/3 -|- C)), 
where Aw is the enthalpy difference and A, 77 and ( are the thermal conductivity and the 




(T/(7)^/^exp {-AttRI/ST) . 
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shear and bulk viscosities, respectively. The droplet formation rate can be expressed as 
(Heiselberg 1995) 

7/2 2 / 3 \ 



with //400 is the quark chemical potential in units of 400 MeV. 

The total number N of droplets formed by nucleation is the integrated rate over the 
volume of the neutron star and the time after the nucleation process starts at the moment 
to, N — J^Anr'^dr I[Ap{r) ,T{t)]dt, where R is the radius of the neutron star. The 
pressure and temperature depend sensitively on the equation of state of the nuclear and 
quark matter. To convert the core of the neutron star into quark matter at least one droplet 
must be formed, i.e., N > 1. Then, with the use of Eq. (3), the condition for the formation 
of at least one droplet requires cr < 24 MeV fm~^AP|(o^^^,io> where Pc and are the core 
values of the pressure and temperature, respectively (Heiselberg 1995). The pressure and 
enthalpy difference depends strongly on the equation of state of nuclear and quark matter 
and hence the droplet formation rate cannot be reliably estimated. 

The effect of subcritical hadron bubbles on an inhomogeneous first order quark-hadron 
phase transition was studied in Shukla et al. (2000). The transition from nuclear matter 
(consisting of neutrons, protons and electrons) to matter containing strangeness was con- 
sidered, within a mean field type model and with the use of Langer nucleation theory, by 
Norsen (2002). An estimate of the time the new phase to appear at various densities and 
times in the cooling history of a protoneutron star was also given. 

On the other hand, as a result of an increase of the central density of the star, due 
for example to accretion or spinning down, a metastable, super-compressed neutron phase 
can appear, with a star consisting of the new phase surrounded by normal neutron matter 
(Grassi 1998). This situation is similar to a gas undergoing a phase transition to its liquid 
state if compressed to a volume Vc at fixed T. If it is slowly compressed, it may stay 
in the vapor state even for V < (similar examples are liquid water freed of dissolved 
air which may be heated above boiling temperature without the formation of vapor and 
cooled below freezing temperature without solidification). The transition from the super- 
compressed neutron phase to the quark phase can take place in a very long time, which in 
some cases can be longer than the estimated fife of the neutron star (Grassi 1998). Therefore 
an increase in the density of the neutron star does not automatically lead to the transition 
to quark matter inside the star. 

The effect of the magnetic field on the quark star structure and on the nucleation process 
of the quark bubbles has been studied in (Ghakrabarty 1995; Ghakrabarty & Sahu 1996; 
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Chakrabarty 1996). In the presence of strong magnetic fields tfie equation of state of strange 
quark matter changes significantly. The strange stars become more compact, the magnetic 
field reducing the mass and radius of the star. The surface energy and the curvature term of 
the quark phase both diverge, with the surface tension diverging logarithmically, while the 
curvature term diverges much faster. Therefore the thermal nucleation of quark bubbles in 
a compact metastable state of neutron matter is completely forbidden in the presence of a 
strong magnetic field. 

These results for the formation of quark bubbles in neutron matter have been obtained 
by using the Cscrnai and Kapusta theory of nucleation (Cscrnai & Kapusta 1992a,b) and 
its extension including the thermal conductivity of dense matter (Venugopalan & Vischer 
1994). In this theory both the hadron and quark materials are considered as substances 
with zero baryonic number, and are treated by using a relativistic formalism. The principal 
result in this approach is the suggestion that the prefactor k is proportional to the transport 
coefficients (viscosity and thermal conductivity) of the neutron matter, and thus when these 
coefficients vanish a quark bubble necessarily does not form. 

However, Ruggeri & Friedman (1996) argued that the energy flow does not vanish in 
absence of any heat conduction or viscous damping. Since the change of energy density e 
in time is given, in the low velocity limit, by the conservation equation de/dt = —V (wv) 
(Csernai & Kapusta 1992a), where w is the enthalpy and v the velocity, this implies that 
the energy flow ~ wv is always present. Therefore, an expression for the prefactor can be 
derived, which does not vanish in the absence of viscosity. The viscous effects cause only 
a small perturbation to the prefactor. The differences between the Csernai-Kapusta (CK) 
(Csernai & Kapusta 1992a) and Ruggieri- Friedman (RF) (Ruggeri & Friedman 1996) results 
are due to the technical differences in the treatment of the pressure gradients. 

A generalized approach, following the CK formalism leading to the prefactor in both 
viscous and non- viscous regimes was developed by Shukla, Mohanty & Gupta (2001). Unlike 
in Csernai & Kapusta (1992a), the linearized relativistic hydrodynamics equations have been 
solved in all regions. 

It is the purpose of the present paper to extend the Zel'dovich nucleation theory for 
the formation of quark droplets in neutron matter, by taking into account both the effects 
of the energy ffow and of thermal conductivity, shear and bulk viscosities. As a result a 
more general expression for the quark matter droplet rate formation can be obtained, which 
is also valid in the limiting case of vanishing conductivity and viscosity coefficients. In the 
nucleation theory of Csernai & Kapusta (1992a), the transition from the neutron to quark 
state is possible only for a normal matter having viscous properties. For zero bulk and shear 
viscosities the transition is impossible. This assumption is, however, to restrictive, and the 
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transition can also take place for a perfect (no viscosity or heat conduction) neutron matter 
fluid. We also derive the expression of the rate formation for this case. 

An alternative point of view of the transformation of neutron matter into strange matter 
can be developed if one assumes that the conversion process is endothermic. In this case the 
strange quark matter formation inside a neutron star is triggered by an external source of 
energy (accretion from the companion star or spin-down). Pulsars are born with an enormous 
angular momentum and rotational energy, which they radiate over a long period of time via 
electromagnetic radiation and electron-positron pair emission. When rotating rapidly, a 
pulsar is centrifugally flattened. With decreasing angular velocity the central density of the 
star is increasing and may attain the critical density necessary for a phase transition. First 
at the center and then in an expanding region the neutron matter will be converted to a 
highly compressible quark matter phase (Glendenning, Pei & Weber 1997). The conversion 
of neutron matter to quark matter alters the moment of inertia of the star, and the epoch over 
which conversion takes place will be signaled in the spin-down characteristics of the pulsar. 
A measurable quantity, the braking index QQ/Q'^, where Q is the angular velocity of the 
star, could be an observational indicator of the slow transition from neutron to quark phase 
(Glendenning, Pei & Weber 1997). By using the formalism developed for the nucleation of 
quark matter we reconsider the idea of phase transition from neutron matter to quark matter 
in rotating compact stellar objects, by assuming that the nucleation and the growth of the 
quark droplets and the formation of a quark core in a neutron star is a result of the transfer 
of the gravitational energy to the quark droplets, due to the spinning down of the rapidly 
rotating neutron star. 

The present paper is organized as follows. The nucleation kinetics of strange matter 
bubbles is considered, within the framework of the Zel'dovich nucleation theory, in Section 
II. The expressions for the rate of formation of quark matter are derived, in both the perfect 
and dissipative neutron fluid case, in Section III. The transition of a neutron star to a quark 
star, due to increase of the central density as a result of the spinning down is considered in 
Section IV. In Section V we discuss and conclude our results. 



2. Nucleation kinetics of strange matter in neutron stars 

We consider that the change from the metastable neutron phase to the stable quark 
phase occurs as the result of fluctuations in a homogeneous medium, formed of neutrons, in 
which small quantities of the quark phase (called bubbles or nuclei) are randomly generated. 
Since the process of creation of an interface is energetically unfavorable, it follows that when 
a quark nucleus is below a certain size, it is unstable and disappears again. Surface effects 
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disfavor the survival of small bubbles below the radius Rc (called critical size-the nuclei of 
this size are called critical nuclei or bubbles) , which is nothing but the value that extremizes 
the thermodynamical work W necessary to create the bubbles. Only nuclei whose size r is 
above the value Rc are stable, and can survive (Landau & Lifshitz 1980). 

The nuclei are assumed to be macroscopic objects containing a large number of particles 
(quarks) . 

For strongly degenerate neutron matter the thermodynamic work W necessary to create 
a quark bubble is (Alcock & Ohnto 1989; Mardor & Svetitsky 1991; Madsen & Olesen 1991; 
Olesen & Madsen 1994; Horvath 1994) 

Air 

W = K - f^n) - (P, - Pn)] yr^ + Anar^ + Sirjr + E^, (4) 

where „ are the pressures of the quark matter and neutron matter (exterior to the bubble, 
assumed to be spherical), a is the surface tension, Ug is the particle-number density of quark 
matter, are the chemical potentials of each phase, 7 is the curvature coefficient and 
Ec = (3/5) Z'^e^ /r is the Coulomb energy of the droplet (Heiselberg, Pethick & Staubo 
1993). Here Ze = {pq — pn) Vd is the excess charge of the droplet with volume Vd compared 
with the surrounding medium, and Pg and p„ are the charge densities of the quark and 
nuclear matter, respectively. Therefore the Coulomb energy of the droplet is given by E^ — 
167r^ {pq — pnf r^/15. Assuming that the quark matter is immersed in a uniform background 
of electros, and that it is in 13 equilibrium, with pd = Ps = Pu + Pe, it follows that the total 
electric charge in the quark phase can be expressed as pq ^ epq (m^/2 — 2pf,pq) /tt^, where 
fJ'u~IJ'd = IJ'q (Heiselberg, Pethick & Staubo 1993). Generally Pg » pn- 

The appropriate form of the curvature energy can be obtained from (Madsen 1993; 
Heiselberg, Pethick & Staubo 1993) 

Ecurv = ^ r [1 + ((^ - /^/^))]~' = ^ [1 + O ^T/pf] , (5) 

where g is the statistical weight and p the chemical potential. Hence the curvature coefficient 
is given by 7 = Zp'^jSi:'^ ~ 18 MeVtm"! (/x/300 MeV)^ (Heiselberg, Pethick & Staubo 1993). 

Requiring W to be an extreme, dW/dr — 0, gives the following equation for the value 
of the critical radius r — Rc- 

IGtt^ 

47r [ug {pg - Pn) - {Pq - Pn)] + Snar + Sn^ = — {pg - pnf r^. (6) 

By neglecting the Coulomb energy, we obtain for the critical radius i?co of the bubble 
the expression 

i^co = ^ (1 + a/1 + &o) , (7) 
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where 
and 



C = -Uq illg - lln) + {Pg " Pn) = -rigA/I + AP > 0, 



6o = 2C|7|/(7'>0. 



(8) 



(9) 



The parameter bo can also be represented as 60 ~ 6 {C/20 MeVfm-2) (7/I8 MeVfrn"^) 
[a/ (75 MeV)^]~^ (Horvath 1994). 

To find the value of the critical radius for a non-negligible electrostatic energy of the 
bubble we use an iterative method, by substituting in the right-hand side of Eq. (6) r by the 
zero'th order approximation given by Eq. (7). Then, in the first order, the critical radius is 
given by 

i?e = ^ (1 + Vl + bj , (10) 

where 



C 



2^+ — {Pg-pnrR, 



2 d4 
■cO 



(11) 



For values of C = 20 MeVfm" 



7 = 18 MeVfm-^ a in the range (75 - 100 MeV)^ 
Pg ^ — 0.4e fm^^ and negligible p„ the critical radius is of the order oi Rc = 2 — 7 fm. 
The Coulomb correction has the effect of increasing the critical radius. The value r = 
corresponds to the limit beyond which large quantities of the quark phase begin to be formed. 
In fact, it is more appropriate to refer not to a limit point r = R^ but to a critical range 

1 /2 

of values of r near that point, with width Sr ~ (T/ATia) ' (Landau & Lifshitz 1980). The 
fluctuational development of nuclei in this range can still, with high probability, throw them 
back into the subcritical range, but nuclei beyond the critical range will inevitably develop 
into the new quark phase. 

The minimum critical work required to form a stable quark bubble is therefore given by 



47rcr^^,,, 



3C2 



where we denoted 



l + 47r (pg - pnf -^G{b, bo) 



F{b) = 2 + 2 (1 + bf^ + 36, 



and 



G(6,6o) 



(i + v^rT6)^ 

2 + 2 (1 + bf'^ + 36 



1 (1 + ^/l + b;) 

5 (i + ^rT6)' 



(12) 

(13) 
(14) 



Within a purely thermodynamic approach, one can put only the problem of calculating 
the probability of occurrence in a medium of fluctuational nuclei of various sizes, the medium 
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being regarded as in equilibrium. Instead of the thermodynamic probabihty of nucleation, 
it is more convenient to use the equihbrium distribution function for nuclei of various sizes 
existing in the medium, denoted fo{r). fodr is the number of nuclei per unit volume of the 
medium with sizes in the range dr. According to the thermodynamic theory of fluctuations 
(Landau & Lifshitz 1980), 

fo{r)^exp[-W{r)/T]. (15) 



Near r — Rc the thermodynamic work can be expressed as 



W{r) 



m 



1 + ^ (p. - Pnf ^,G{b) 



UTraC^Vm - 327r2 (p^ - p^y (l + ^/TTb} 

3^3 



{r-Rcf. (16) 



Therefore for the equilibrium distribution function we find the expression 
'UnaC^y/T+b - 3277^ {pg - p^f (l + VT+b)^ 



fo{r) = fo (Rc) exp 



3C^T 



T - R,y 



where /o {Rc) = Cq x exp (— Cq =constant. 



(17) 



In order to estimate the coefficient Co of the exponential in /o (-Rc) we follow Landau & 
Lifshitz (1980), and assume Co = UgUnRl, with Un and Ug the particle number densities in 
the neutron and quark phases, respectively. Thus we obtain 



fo (Rc) = rigTinRc exp 



l + 47r(p,-pj^— G(6,M 



(18) 



Let /(t, r) be the kinetic size distribution function of the nuclei. The elementary process 
which changes the size of a nucleus is the attachment to it, or the loss by it, of a quark 
droplet, and this is to be regarded as a small change, since the nuclei arc considered to be 
macroscopic objects. Therefore the growth of the nuclei is described by a Fokker-Planck 
equation (Landau & Lifshitz 1980) 

where j = —B^ + Af is the flux in the size space. B is the nuclear size diffusion coeffi- 
cient and A is connected with i? by a relationship, which follows from the fact that for an 
equilibrium distribution j = 0. Therefore we find A = —BW'{r)/T. 
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In the case of a continuous stationary phase-transition process we have j ^constant. 
The constant flux is just the number of nuclei passing through the critical range per unit 
time per unit volume of the medium, i.e. it defines the rate of the process. With the use of 
the condition of the constant flux we obtain —BJq-^ = j giving 

f f dr 

f = -J / ^ + const. (20) 

Jo J -D/O 

The constant in this equation and j are found from the boundary conditions for small 
and large r. The fluctuation probability increases rapidly with decreasing size and small 
nuclei have a high probability of occurrence. This is expressed by the boundary condition 
///o— >lasr— >0. The boundary condition for large r can be estabhshed by noting that 
above the critical range the function /o increases without limit, whereas the true distribution 
function /(r) remains finite. This situation is expressed by imposing the boundary condition 
///o = for r — > oo. The solution which satisfies the above conditions is (Landau & Lifshitz 
1980) 

/ ,r^.i.r* pi) 



fo Jr Bfo j Jq Bfo 

In these equations the integrand has a sharp maximum at r = Rc. By extending the 
integration with respect to r from — oo to +oo, one obtain for the number of viable quark 
nuclei formed in stationary conditions per unit time and per unit volume the expression: 



Above the critical range, the distribution function is constant: having reached that 
point, the nucleus becomes steadily larger, with practically no change in the reverse direction. 
Accordingly we can neglect the term containing the derivative df / dr in the fiux, leading to 
j — Af. Prom the significance of the fiux it follows that the coefficient A acts as a velocity 
in size space, A — {dr/dt)^^^^ (Landau & Lifshitz 1980). Therefore we find for B: 

T f df 



B 



W'{r) \dt 

^ ..3) 



UnaC^VTTb - 32%^ {pg - pnf (l + VTTS) (r - Re) 
The rate of growth of the bubble radius near the critical radius is given by (^) = 

° ° \dt J macro 

k{r — Rc), where k is the dynamical prefactor (Langer & Turski 1973; Csernai & Kapusta 
1992a). 
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The prefactor k has been evaluated, by solving the equations of relativistic fluid dynam- 
ics in all regions, in Shukla, Mohanty & Gupta (2001). The result (also taking into account 
heat conduction) is 

where Cg is a constant (the velocity of the sound in the medium around the saddle con- 
figuration) and A„, r]n and ^„ are the thermal conductivity and shear and bulk viscosity 
coefficients of the neutron matter, respectively. The first term in the above equation is the 
same as obtained by Ruggeri & Friedman (1996), corresponding to the case of non- viscous 
matter. The second term is similar to the results obtained by Csernai & Kapusta (1992a) 
and Venugopalan & Vischer (1994), except with a minor difference , i.e. instead of 4 there 
is a factor in the numerator. 



(24) 



3. Neutron matter to queirk matter transition rates 



With the use of Eq. (24) for the prefactor it follows that the rate of formation of quark 
bubbles in neutron matter is given by 



3 = 



Y 12aC^VT+b - 327r (p^ - p„)' (l + VT+b) ' 
{XnT + 2 {lr]n + Q) 



Rl {^wf 



+ 



clRl {AwY 



'^F{b) l + 47ripg-pnT^G{b,bo 



exp 



2 ^2 



T 



(25) 



To obtain the number ^ of the net strange quark matter bubbles, the rate j must be 
multiplied by the time interval available for prompt nuclcation. At and by the volume Vq, 
where the nucleation can take place in the dense core (Horvath, Benvenuto & Vucetich 1992). 
We impose that at least one quark bubble appears (which would suffice to convert the whole 
neutron star). Therefore, with the use of Eq. (25), one obtains the following general condition 
for a quark bubble formation in a neutron star core: 



e = jAtAVo > 1. 



(26) 



The quark bubble formation rate essentially depends on the enthalpy difference in the 
two phases. In the density range of interest, all the neutron matter EOS can be very 
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well parameterized as polytropes, P„ = Kri^^ (Lugones, Benvenuto & Vucetich 1994). The 
corresponding energy density in the neutron phase is = UBrrin + :::^Kn1^ , leading to 
Wn — Sn + Pn = nsTUn + ;^i^n^. We also assume that the formed quark bubble consist of 
u and d quarks in the ratio 1 : 2; only later weak interactions may change the composition 
to an energetically more favorable state. The quarks chemical potentials are related by 
fid = 2^^^ flu, and, assuming chemical equilibrium across the phase boundary, we also have 
l^n = l^u + '^l^d = (1 + 2^/^) flu (Olesen & Madsen 1994). Then the pressure in the quark 

. 4 _|_ ..4 

phase is (assuming a simple bag model) Pq = ^ 2'* — B and Eg — 3Pg + B (Cheng, Dai & 



Lu 1998), giving Wq = ^"^^'^ — ZB. Hence, the enthalpy difference in the two phases can be 
approximated by 

A A 

7 



Aw = 



4 I 4 
7r2 



3B — nBrrin 



7 



(27) 



Assuming that the energy flow is provided by the viscous effects only, one obtains for 
the pref actor the expression (Venugopalan & Vischer 1994) 



2a 



A„r + 2 I -r?„ + e„ 



(28) 



In the limit of zero baryon number, A„ — > and we obtain the result of Csernai & 
Kapusta (1992a). If the matter is baryon-rich, but viscous damping is negligible, ?7„, ^„ — > 0, 
and we obtain the results of Langer & Turski (1973) and Kawasaki (1975). 

Therefore the condition of the formation of a quark bubble is given by 
1 n^rin / ZC^Ta 



TT (A«;) R, y l2C^,/lTh - 327r (p, - pnf (l + ^/TTiy 
F{b) l + 47r(p. 



A„T + 2( -ry^ + Cn 



exp 



Pnf ^G{b, bo) 



X 



AtAVn > 1. 



(29) 



The transport properties of dense matter have been intensively investigated in both 
high energy physics and astrophysical frameworks (Flowers & Itoh 1976, 1979, 1981; Hakim 
& Mornas 1993). The shear viscosity rj and thermal conductivity A for neutron mat- 
ter have been derived by Danielewitz (1984) from the Uhlenbeck-Uehling equation (see 
also Sawyer (1989) and Cutler, Lindblom & Splinter (1990)). They are given by 77 = 
{16n^)~' {p%/mlai) T'"^ and A = 5 (pf^/m^aa) r-796, where PiT is the Fermi momentum 
and CTi, a'2 are some quantities related to particle-particle cross sections, estimated in Daniele- 
witz (1984). The contribution of the bulk viscosity and thermal conductivity are negligible, 
and the main contribution is from the shear viscosity, which is typically of the order of ~ 50 
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MeV/fm^ (Heiselberg 1995). But in the case of an ideal neutron gas, with zero viscosity, 
there will be no bubble growth, and in this case the transition from neutron matter to quark 
matter in astrophysical objects cannot take place. 

Assuming that the energy flow does not vanish in the absence of any heat conduction 
or viscous damping, and considering that the viscous effects are small and can be neglected, 
the condition of the formation of a quark bubble inside the dense core of a neutron star 
becomes: 



1 nqUn 



v^TT Aw Y uC^y^T+b - 32n {pg - p^f cr^ (l + y/lTb) 



X 



exp 



3C^T 



l + 47r(p,-pO'^G(6,6o) 



AtVn > 1. 



(30) 



For r > Rc the time growth of the radius of the quark droplets can be described by the 
equation (Csernai & Kapusta 1992a; Shukla, Mohanty & Gupta 2001) 



with the general solution 
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(31) 



(32) 



where we have used the initial condition r(0) = 2i?c. Therefore for a neutron star with 
radius R, the time tconv required for the conversion of the whole star to a quark star can be 
obtained from 
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(33) 



Since R» Rc, with a very good approximation we find 

2 



tc 



1 R 



2k \R. 



(34) 



Neglecting the viscous effects, and by assuming again that the neutron matter is de- 
scribed by a polytropic equation of state, while the quark phase obeys the simple bag model 
equation of state, one obtains for the conversion time 



tr 



\^8aRcW, 
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(35) 
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The volume V{t) of the neutron matter converted to quark matter after a time interval 
t < tconv is given, in this simple model, by 



3 (Aw) ' 



(36) 



If the transition from neutron to quark phase is driven by viscous processes only, then 



f ^ R 



\ 2 



2a [A„r + 2(|77„ + e„)] 
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(37) 
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4. Endothermic nucleation of quark matter bubbles in the core of rotating 

neutron stars 

In the previous Sections we have considered the process of nucleation of quark bubbles 
in neutron star cores, and we have derived the conditions for quark matter formation in an 
astrophysical context. In the present Section we shall apply the results previously obtained 
to analyze the possibility of the transition from neutron to quark phase during the spinning 
down of a rapidly rotating neutron star. In this case the rotational energy of the star can 
be used to trigger the phase transition. 

The minimum energy density required to create a quark bubble is 



C 

Ua — - — — (0) 



247r , ,2 o-^ 

1 + — (p,-Pn)'^G(&) 



(39) 



An estimation of Eq. (39) for C ^ 20 MeVfm^^ 7 f« 18 MeVfm-\ a ^ (75 MeV)^ 
Pq Ri — 0.4e fm~^ and Pn ~ gives Ug ^ 3.76 x 10^^ erg cm~^. 

We consider that the evolution of the rotating neutron star can be approximated by a 
sequence of MacLaurin spheroids (Chandrasekhar 1986; Cheng et al. 1992). The mass of the 
star will be denoted by M, and the major and minor axis of the star by a and c (equatorial 



and polar radius). Then the eccentricity of the star is defined according toe^yl — (^) . 
Let p ^constant denote the density of the star, corresponding to a given value of the angular 
velocity Q and of the equatorial radius a. 
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Hence the basic equations describing the mass, hydrostatic equihbrium and gravitational 
energy Eg^ of the rotating star are (Chandrasekhar 1986) 

3 GM^ arcsin e 



47r 



where 



M = —a'p (1 - , = 2nGpf{e), Egr - 

O O (X G 



^, , il-eY\^ o 2x . 3(1-6^) 
/(e) = — (3 - 2e^) arcsme — - — -. 



(40) 



(41) 



Assuming that the total mass of the star is a constant, one obtain the following relation 
between the variation Aa, AQ, Ap and Ae of the equatorial radius, angular velocity,density 
and eccentricity, respectively: 



3 Aa 1 Ap _ e Ae 
a AQ ^ p AQ " 1 - e2 AQ ' 



(42) 



The variation of the eccentricity with respect to the angular velocity is given by 

Ae _Q-7rG'/(e)^ 



where 



An TrGpg{e) ' 
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With the use of the above equations the variation of the gravitational energy of the 
neutron star, corresponding to a simultaneous change in equatorial radius, angular velocity 
and density, is given by 



AEgr (a, p, n) 



3GM^ 
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(45) 

(46) 
(47) 



The gravitational energy density released as the result of the slowing down of the neutron 
star can trigger the endothermic phase transition at the center of the star. Inside a sphere 
of radius Rq the gravitational energy density is 

9 GM^ 



Ugr (a, p, Q) 
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(48) 
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This energy should be greater or equal to the minimum energy necessary for the for- 
mation of a quark bubble, given by Eq. (39). Therefore the gravitational energy converts to 
the quark phase the neutron matter inside a sphere of radius 



Rq{a,p,n) = 2.88 



M 
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(49) 



In order the whole neutron star be converted to a quark star the condition Rq {p, Q) = 
a (p, Q) must hold. On the other hand, for typical neutron star densities of the order of 
10^"^ — 10^^ g/cm"^, the term A (e) Q/nGp is much smaller than B (e) (Ap/AQ) /p, 



n „ , , 1 Art 



(50) 



For small eccentricities e the function B{e) can be approximated by 1/3. 

Therefore we obtain the following condition for the relative change in the central density, 
necessary to convert a neutron star to a strange star, due to the endothermic gravitational 
energy transfer in the whole volume of the star: 



Ap 
P 



C 



0.125 X 



2 + 2(1 + 6)^/^ + 36 



M 



-2/3 



(51) 



For a neutron star of around two solar masses, M 
MeV and b ^ 6, we obtain ^ ^ 1.85. 



2Mq and by assuming C ^ 20 



Of course this result is strongly dependent on the precise numerical values of C and b, 
which are generally poorly known. C is also a neutron matter equation of state dependent 
parameter. For example, assuming for C a value of 10 MeV will lead to ^ ^ 0.92. An 
increase in the mass M of the star will also decrease the value of the relative change in the 
density of the star necessary to convert the neutron star. 

At present, we have very little observational knowledge about how fast newborn neutron 
star can rotate. The most well-known young pulsar is that the Crab pulsar was born with a 
rotational period of about 20 ms (Manchester & Taylor 1977). If so, during its hfetime the 
pulsar its central density only increases about Apd Pc ~ 0.001 (Ma &; Xie 1996). Currently 
the fastest rotating pulsar known is PSR1937+214, which has a period of 1.55ms or Q 
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4, 000s~^. But it has weak magnetic field and is suggested to be spun-up by accretion (Alpar 
et al. 1982). However, some rapidly spinning millisecond pulsars are suggested to be born 
by accretion induced collapse from white dwarfs (Arons 1983). Therefore it can not ruled 
out that some pulsars can be born with milliseconds and strong magnetic field. 

The theoretical investigation of rotating general relativistic objects performed by Cook, 
Shapiro & Teukolsky (1994) shows that for some realistic equations of state of neutron matter 
this variation of the central density can be achieved during the complete spin-down of the 
star (see Table I). 



EOS 


Mstot(in Mq units) 


Rstat (km) 


An (s-1) 


Pc 


A 


1.6551 


8.368 


10011 


1.26 


Ai; 


2.1335 


9.111 


10587 


0.97 


FPS 


1.7995 


9.281 


8874.9 


0.91 


L 


2.7002 


13.7 


6482.9 


1.244 


M 


1.8045 


11.6 


4437 


4.189 



Table 1: Variation of the central density of the rapidly rotating neutron star for different equa- 
tions of state (Cook, Shapiro & Teukolsky 1994): A (Pandharipande 1971), AU (Wiringa, 
Fiks & Fabrocini 1988), FPS (Lorentz, Ravenhall & Pethick 1993), L (Pandharipande Pines 
& Smith)) and M (Pandharipande & Smith 1975). 

The central density increase due to the spin-down can be easily realized in a short 
time interval for a special class of stellar type objects, called magnetars. Magnetars are 
compact objects with super strong magnetic fields of the order B fa 10^^ Gauss or even 
higher (Duncan & Thompson 1992; Paczynski 1992; Thompson & Duncan 1995, 1996). It 
is now believed that soft gamma repeaters (SCR's) - a small class (4 confirmed and one 
candidate) of high energy transient discovered through their emission of bright X-ray/7-ray 
bursts, which repeat on timescales of seconds to years - are magnetars. There is evidence that 
the giant SCR flares involve the cooling of a confined e^-photon plasma in an ultrastrong 
magnetic field (Lyubarsky, Eichler & Thompson 2002). Some authors (Cheng & Dai 1998; 
Usov 2001) have suggested that SCR's are strange stars, motivated in part by the super- 
Eddington luminosities of their giant flares. The magnetars differ from the canonical pulsars 
(with low magnetic fields of the order B 10^^ — 10^^ Causs) in the sense that they spin 
down much more rapidly. 

If we assume that the spin-down of the magnetar is completely determined by the torque 
of its relativistic wind emission, generated via the magnetic dipole radiation, then the time 
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variation of the angular velocity Q of the star is given by Shapiro & Teukolsky (1983) 



(52) 



where / is the moment of inertia of the star and /j, = R^B is the magnetic dipole moment. 
When the star is born with a spin period much shorter that the observed one, the age of 
the star is the spin-down age Tsd — P/ P, where P and P are the present spin period and its 
time derivative. If the spin-down is entirely due to the magnetic dipole radiation we obtain 

_n_ 3 ic' 

By adopting for the moment of inertia and the radius of the magnetar the typical values 
/ = 10"^^ gcm^ and R = 10^ cm, and by assuming that the star was born with a magnetic 
field of the order of S ~ 10^^ Gauss and with an angular velocity of Qq = 6000 rads~^, the 
spin down age is given by Tgd ~ 1125 s. 

From Eq. (52) it follows that the decay law of the angular velocity fl of the star is of 
the form 

!2 = , "° (54) 

where is the initial angular velocity of the star. 

Hence, in a time interval of around t ~ 10^ s and for a magnetic field B = 10^^ G the 
angular velocity of the star decreases from Qq = 6000 rads^^ to ^2 = 142 rads^^. For an 
initial angular velocity of the order of f2 = 8000 rads^^ and for a magnetic field of the order 
of -B ~ 10^^ G, Tsd ~ 5 s. In a time interval of around t ~ 600 s the angular velocity of the 
star decreases to ^2 fti 579 rads~^. 

Therefore in a magnetar density changes, due to the spin-down, take place in a short 
interval of time. Hence they can provide enough energy to trigger spontaneous nucleation 
to the quark phase, once the central density of the star goes above the critical density. 

The nucleation of quark droplets in the presence of strong magnetic fields essentially 
depends on the strength of the magnetic field. The energy of a quantum particle changes 
significantly if the magnetic field is of the order or exceeds the critical value Be = m'^c^/eh, 
where m is the mass of the particle. For electrons the value of the critical field is i?c^^ = 
4.4 X 10^^ G. In order to calculate the value of the critical magnetic field for quarks we have 
to use the estimations for the quark masses. However, there are many uncertainties in the 
numerical values of the quark masses, which play an essential role in the calculation of the 
critical magnetic field and which are generally poorly known. In the case of free quarks. 



-20- 



the current quark mass is considered to be in the range of ruq ~ 5 — 10 MeV. Then the 
corresponding critical magnetic field is Bc'^^ ~ (1 — 2) x 10^ x B^'^^ (4.4 — 8.8) x 10^^ G. 
For the current quark mass, this is the typical strength of the magnetic field at which the 
cyclotron lines begin to occur. In this limit the cyclotron quantum is of the order or greater 
than its rest energy. This is equivalent to the requirement that the de Broglie wavelength is of 
the order or greater than the Larmor radius of the particle in the magnetic field (Chakrabarty 
1995). However, from the physical point of view a better description of the quark mass can 
be obtained by using instead the current quark mass the effective or constitute quark mass, 
which takes into account the effects of the strong interactions, and which is of the order 
rriq ~ 100 — 300 MeV (Cheng, Dai Sz Lu 1998). The critical magnetic field corresponding to 
the effective quark mass could be as high as sl^^ ~ 10^^ - 10^^ G. 

The study of the behavior of the surface and curvature energy terms in strong mag- 
netic fields have shown that both quantities diverge for B > B'f^, with the curvature term 
diverging much faster. Consequently, in the presence of strong magnetic fields the rate of 
stable quark droplet formation per unit volume / ~ T'^cxp [^W/T) tends to zero, 7 — > 0, 
and therefore there cannot be any thermal nucleation of quark droplets in neutron star cores 
(Chakrabarty 1995, 1996). 

The value of B inside a magnetar is not known, and it is not certain if the interior 
magnetic field is much stronger than the surface field. If in the interior of the magnetar the 
value of the magnetic field can exceed the critical value B^'^ ~ 5 x 10^^ G, which represents 
the most conservative estimation of the critical magnetic field, corresponding to the free 
quark mass m = 5 MeV, then the rate of droplet formation can be considerably reduced and 
the transition from neutron matter to quark matter could not take place. 

On the other hand, by taking into account that the physical mass of the quark is 
the effective mass, the corresponding critical magnetic field could be of the order B^/^ ~ 
10^^ — 10^^ G, value which could be higher than the magnitude of the magnetic field inside 
the magnetars. For magnetic fields smaller than B^c \ the endothermic nucleation process 
can take place in a short time interval during the spin-down of the magnetar. 

5. Discussions and final remarks 

In the present paper we have considered a kinetic theory of strange matter nucleation 
in neutron stars. We derived general necessary conditions for neutron to quark matter 
conversion. There have been several assumptions made in the derivation of the results, 
mostly related to the form of the growth speed of the bubble. The result is vahd only if 
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non-linear effects can be ignored, and tfie linearized hydrodynamic equations are applicable. 
Furthermore, we have assumed that the radii of the bubbles are larger than the correlation 
length and that heating due to dissipation is slow, causing the temperature to vary slowly 
across the bubble wall. We also considered that the phase transition is strongly first order, 
releasing considerably heat. 

The growth of the critical size quark bubbles nucleated in the first order transition in the 
neutron star core is governed by the dynamical prefactor k. There are several expressions for 
the prefactor obtained in the physical literature. In the Cscrnai-Kapusta approach (Csernai 
& Kapusta 1992a), the prefactor essentially depends on the dissipativc properties (shear, 
bulk viscosity and thermal conductivity) of neutron matter. Therefore, in a non-viscous 
medium with in = Vn = = 0, the growth of the quark bubbles cannot take place. The 
main physical assumption is that the energy flow is provided by the viscous effects only. 

On the contrary, one can argue that the energy flow does not vanish in absence of 
heat conduction or viscous damping (Ruggeri & Friedman 1996; Shukla, Mohanty & Gupta 
2001). In this case, and in the hmit of zero viscosity, the prefactor depends only on two 
scale parameters, the correlation length ^ and the critical radius of the quark bubble Rc- 
For a viscous medium the prefactor can simply be written as the sum of the viscous and 
non-viscous term, with the viscous term not affecting the growth process significantly. The 
previous (viscous) results can be rcobtaincd by using some assumptions for the velocity of 
sound in the medium around the saddle configuration. 

The condition for the formation of at least quark bubble in the neutron matter as a 
result of the thermal fluctuation is given by Eq. (29). The condition includes both the 
cases of viscous and ideal neutron matter. It is extremely sensitive with respect to the 
numerical values of the microscopic parameters characterizing the quark bubble, like the 
surface tension a, the bag constant B or the curvature coefficient 7. In discussing the 
astrophysical implications of the expression (29) we shall consider separately the two limiting 
cases, corresponding to viscous and non- viscous neutron matter, respectively. 

In the case of ideal neutron matter, the condition of the formation of a quark bubble is 
given by Eq. (30). This equation determines at which temperature the phase transition takes 
place. For example, in the case of the Walecka mean field equation of state (Walecka 1975), 
with K — 15758 and 7 = 4.95 (Lugones, Benvenuto & Vucetich 1994), the temperature Tc 
necessary for the formation of a stable quark bubble in a second {At = 1 s), in a volume 
AT^o with radius ARq — 1 km {AVq — Att/S km^), and at the center of the neutron star, is 
around Tc = 12 — 13 MeV, for a = (85 MeV)^. Here we have also assumed the standard 
values for the quark matter chemical potential and for the bag constant, yU„ = = 280 
MeV and B = 60 MeVfm"^ (Cheng, Dai & Lu 1998). Since the temperature of a newborn 
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neutron star is around 10 MeV (Olesen & Madsen 1994), this value of the surface tension 
seems to rule out, for the given equation of state, the possibility of formation of quark nuclei 
in the neutron matter. However, a small variation in the value of a, ex = (75 MeV)^ reduces 
the temperature for quark nuclei formation to = 4 — 6 MeV, a value which docs not 
exclude the possibility that quark nuclei formation can be initiated in the early stages of the 
neutron star evolution. Smaller values of a, of the order of cr = (50 MeV)^, can lower even 
more the phase transition temperature, making it possible even for cold neutron stars. The 
transition temperature is relatively insensitive to the details of the equation of state of the 
dense neutron matter. 

If the temperature at which the phase transition is initiated is strongly dependent by the 
microscopic model adopted for describing the neutron matter properties, the total transition 
time of the neutron star to a quark star, given by Eq. (35), is much less dependent on 
the details of the equation of state or the exact numerical values of the surface tension or 
curvature of the bubble. Generally, for a neutron star with radius i? = 12 — 15 km, the 
conversion time is very long, of the order of 10^ years for a ~ (75 MeV)^. For a ~ (50 
MeV)^ or even lower, the conversion time could have values of the order of 10^ years, which 
are still of the same order of magnitude as the lifetime of a neutron star (10^ years). Hence, 
in the present approach, which is limited to the consideration of linear effects only in the 
hydro dynamic description of the first order phase transitions for an ideal neutron fiuid, 
the growth of the quark nuclei is very slow. Therefore the results of the present analysis 
suggest that even quark nuclei could appear inside the neutron star at a very early stage, 
the transition to a pure quark phase ends only in the final stage of its evolution. 

If the neutron matter-quark matter phase transition is mainly driven by viscous pro- 
cesses, like in the Csernai-Kapusta scenario, then the time necessary for a neutron star to 
convert to a quark star follows from Eqs. (28) and (34). In order to numerically estimate the 
conversion time, we need an estimation of viscosity and conduction coefficients of nuclear 
matter. We shall also consider two limiting cases for the phase transition. Below the critical 
density p = Pc ~ 10^^ g/cm^, cold neutron matter consists of two distinct phases. At low 
density, neutron star matter consists of nuclei which form a solid lattice for temperatures 
below the melting temperature and a sea of relativistic electrons. When the density reaches 
p 10^^ g/cm^, neutrons begin to "drip" from the neutron-rich nuclei and then, in addi- 
tion to the relativistic electrons there is a sea of non-relativistic neutrons. At the critical 
density the neutron-rich nuclei dissolve, leaving seas of degenerate neutrons, protons and 
electrons. The neutrons in cold neutron star matter are quite likely superfiuid, and therefore 
the transport coefficients are dominated by electrons, with the motion of electron deter- 
mined by electron-electron and electron-proton scattering (Flowers & Itoh 1979)As a result 
the viscosity is much lower, but the thermal conductivity is only slightly smaller as compared 
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to the case in which both neutrons and protons form normal fluids (Flowers & Itoh 1979). 
Hence, we can neglect the viscosity coefficients with respect to the thermal conductivity and 
approximate A„ by A„ ~ Xopi^/Tg, with Aq ~ 10^^ (Flowers & Itoh 1979). A slightly different 
expression for the leading term in thermal conductivity has been proposed in Danielewitz 
(1984), A„T ^ Ao (n/no)^^^, with Aq = 0.15 fm~^ and no = 0.145 fm~^. Hence, by using this 
last expression, we obtain for tconv 

f ^ P (^^)' P2 

4(7 Ao [n/no) ' 



The transition time is independent of the temperature of the neutron matter. It depends 
only on the microscopic properties of the quark bubble (surface tension and curvature coef- 
ficient), density and the equation of state of the matter inside the star. For Hu = l^d = 280 
MeV and 5 = 60 MeVfm"^ and for a (75 MeV)"^ Eq. (55) gives, for a neutron star with 
radius i? = 12 km and for several equations of state values of the conversion time of the 
order tconv ~ lO"*^^ years. 

In the case of young and hot neutron stars, the main dissipative mechanism which could 
drive the phase transition is the shear viscosity of the neutron matter. The temperature de- 
pendence of the shear viscosity coefficient can be approximated by 77^ ~ {rjo/T'^) (n/no)^, 
with 770 = 1700 MeV^/fm^ (Danielewitz 1984). Consequently, by neglecting the heat con- 
duction and the bulk viscosity, the conversion time becomes 

tconv ~ ■^-^c'^ ^—R^ ~ —Rc — — 7 — 72"^^^- (^6) 

cr^n <J^ am (n/no) 



In small temperatures, the shear viscosity becomes very large. Hence, for a shear vis- 
cosity driven phase transition the conversion time could be very small. 

In Fig. 1 we have represented tconv given by Eq. (56) as a function of temperature, for 
different equations of state of neutron matter. We have considered four equations of state, 
namely, the equation of state (EOS) of the free neutron gas (Shapiro & Teukolsky 1983), the 
Bethe- Johnson EOS (Shapiro & Teukolsky 1983), the Lattimer-Ravenhall EOS (Lattimer 
& Ravenhall 1978) and the Walecka EOS (Walecka 1975). In all cases we calculated the 
conversion time for a neutron star with radius R — 10 km. 

As one can see from the figure, in this case the conversion time is much shorter as 
compared to the ideal neutron matter case or thermal conduction driven phase transition. 
tconv is rapidly decreasing with the temperature. For a neutron star with rapidly cooling 
processes, e.g. direct URCA process, kaon condensation etc., the stellar temperature can 
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Fig. 1. — Variation of the conversion time tconv for a neutron star with radius R — 10 km, for 
different equations of state of the neutron matter: free neutrons EOS (sohd curve), Bethe- 
Johnson EOS (dotted curve), Lattimer-Ravenhall EOS (short-dashed curve) and Walecka 
EOS (long-dashed curve). 
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decrease to less than 10~^MeV in a time scale of 10~^ yr (Tsuruta 1998). According to 
Fig. 1, the conversion can take place in a time scale less than 10~^ yr if the spin-down 
time scale is less than this time scale. This can be case if this neutron star is a magnetar. 
We should point out that the interior of magnetars are not quantum liquid, i.e. superfiuid 
or superconducting, because the strong magnetic field can destroy the nucleon cooper pairs 
(Chau, Cheng & Ding 1992). 

On the other hand, if the neutrons and protons are superfiuid, then the transport proper- 
ties are determined completely by electron-electron scattering. In fact in this case the viscos- 
ity of the nuclear fluid is given by r^"^ = (157r^Q;^/2pg) (2pe/^/t) [5/2 -|- 3 (me/pe)^ + (^^e/Pe)^] 

1 /2 

(Flowers & Itoh 1976), where a is the fine structure constant, kft = {iapeeg/TT) ' is the 
Thomas-Fermi wave vector, and rUe, Pe and eg are the electron mass, momentum and energy, 
respectively. 

We assume that the electron gas inside the neutron star is extreme rclativistic and 
strongly degenerate. Then it follows that Pe ^ Pf = (3ne/87r)^^^ (Shapiro & Teukolsky 
1983), with pf the Fermi momentum and Ue the electron number density. Since Ug can 
be related to the total particle number n by means of the relation Ue — YeU, with Yg the 
mean number of electrons per baryon, it follows that Pe ~ (3y^no/87r)^''^ {n/rioY^^. Hence 
we obtain for the shear viscosity coefficient of the nuclear matter 
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The variation of the conversion time, as a function of temperature, is represented, for a 
superfiuid neutron and proton core of the star with = 0.1, in Fig. 2. 

In this case, since the viscosity of the superfiuid nuclear matter is much smaller, the 
conversion time is longer than for a nuclear fiuid dominated by the neutron viscosity. How- 
ever, in this case the phase transition from neutron matter to quark matter can also take 
place at very low temperatures. Consequently, a very short tconv, of the order of hours or a 
few days, is also allowed in the framework of this model. 

It is widely accepted that at asymptotic densities the ground state of QCD with nis — 
is the color-fiavor locked (CFL) phase (Alford et al. 1998; Rajagopal Sz Wilczek 2001). In this 
phase, color gauge symmetry is completely broken, as are both chiral symmetry and baryon 
number. The effective coupling is weak and the low-energy properties can be determined by 
using methods from the theory of superconductivity. For large gaps (A ~ lOOMeV), the 
CFL phase is rigorously electrically neutral, despite the unequal quark masses, and even in 
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log[T](MeV) 

Fig. 2. — Variation of the conversion time tconv for a neutron star with radius R— 10 km to 
a quark star, in the case of the presence of a superfluid core, for different equations of state 
of the neutron matter: free neutrons EOS (sohd curve), Bcthc- Johnson EOS (dotted curve), 
Lattimer-RavenhaU EOS (short-dashed curve) and Walecka EOS (long-dashed curve). 
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the presence of electron chemical potential. The transition from nuclear matter to quark 
matter via bubble nucleation may be greatly simplified if the transition occurs directly to 
quark matter in the CFL phase (Rajagopal & Wilczck 2001). For a given baryonic chemical 
potential /i, electrically neutral nuclear matter and electrically neutral quark matter have 
different values of the electronic chemical potential fie- Since fig must be continuous across 
any interface, a mixed phase region is formed, within which positively charged nuclear matter 
and and negatively charged quark matter with the same /ig coexist at any given radius. The 
growth of the quark droplets is at the expense of the nuclear matter. However, if the quark 
matter is in the CFL phase, an interface between bulk nuclear matter with non-zero electron 
number ^ and CFL quark matter with — may be stable, as long as fx^ satisfies the 
condition \ml/Aji — Sfi] < A/\/2, with Sfi the variation of the chemical potential (Rajagopal 
& Wilczek 2001). Therefore, the existence of the CFL phase supports the idea that in 
neutron stars quark matter and baryonic matter may be continuously connected. Then the 
growth of the interface between the neutron and CFL quark phase could be the result of 
exterior energy absorption, via an cndothermic mechanism. 

In summary, by using the kinetic nucleation theory we have considered the possibility 
that the neutron matter-quark matter phase transition could be triggered by the gravitational 
energy released during the spin-down of a pulsar. This process could take place in a very 
short time interval, of the order of few tens seconds, in a class of dense compact objects, 
with very high magnetic fields, called magnetars. However, the percentage of neutron star 
matter that can be converted into quark matter depends on the initial period, the neutron 
star mass and on the equation of state of the neutron and quark phases. An important and 
open issue is how could we determine if such a process is taking place or took already place 
in some pulsars or magnetars, and what are the observational signatures of the endothermic 
neutron matter-quark matter transition. In particular, it would be interesting to know if 
such a process could be a possible gamma-ray burst mechanism, as suggested by Berezhiani 
et al. (2003). All of these issues will be considered in a subsequent paper (Tang, Cheng and 
Harko 2004). 

This work is supported by a RGC grant of the Hong Kong government of the SAR of 
China. The authors would like to thank to the anonymous referee for comments that helped 
to improve the manuscript. 
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